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1.0 TECHNICAL SUMMARY 

The lake acidification in Northern Ontario has been investigated 
using Landsat TM to sense lake volume reflectance and also to provide 
important vegetation and terrain characteristics. The purpose of this 
project was to determine the ability of Landsat to assess water qual- 
ity characteristics associated with lake acidification. Our basic 
hypothesis is that seasonal and multi-year changes in lake optical 
transparency are indicative of reaction to acidic deposition. Results 
from this study demonstrate that a remote sensor can discriminate lake 
transparency based upon measured reflectance. In many acid sensitive 
lakes, optical transparency is controlled by the amount of dissolved 
organic carbon (DOC) present. DOC is a strong absorbing non- 
scattering material which has the greatest impact at short visible 
wavelengths including TM band one. Acid sensitive lakes have high 
concentrations of aluminum, which have been mobilized by acidic com- 
ponents contained in the runoff. Aluminum complexing with DOC is 
considered to be the primary mechanism to account for increased lake 
transparency. 

When eco-physical properties developed from vegetation, soil/ 
bedrock, sulfate deposition, and topographic relief characteristics 
were stratified across the study regions, it was determined that these 
regions could be described as ten separate environments based upon a 
simple acid sensitivity index model. This classification of the 
environment predicts location of regions containing acid sensitive 
lakes. The spatial co-occurrence of acid sensitive eco-physical 
parameters showed that acidification of a lake is driven mostly by 
local geology and soil conditions and less by the rate of sulfate 
deposition. Geologies which are weather resistant containing quartz 
rich sandstones and other quartz rock with bare or shallow sandy soils 
are most susceptible to regional acid deposition. These geologies 
produce naturally very low buffered acid sensitive lakes, contain very 
low amounts of DOC, and tend to have lower values of pH. 
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This study involved gathering an extensive amount of supporting 
data from 1986 and 1987. During August 1986, data were gathered from 
several sites representative of the range of ecosystems found in 
Northern Ontario. These data include limnological parameters, subsur- 
face spectral irradiance, subsurface beam attenuation, airborne 
radiometry, and Landsat TM coverage. Based on these data, lake 
reflectance was modelled in terms of DOC and chlorophyl 1-a pigment 
concentrations. It was demonstrated that acid lakes having abnormally 
small amounts of DOC show greater reflectance than lakes with normal 
pH and DOC values. Significant correlation was found between in-situ 
and above surface lake volume reflectances. The model -predicted chan- 
ges in TM band one signal response were consistent with observed 
values. 

A second data set was gathered during May and June of 1987 on 
eight lakes to observe possible seasonal changes in subsurface and 
Landsat TM reflectance measurements. It was expected that spring 
runoff would produce decreases in DOC concentration and an increase in 
reflectance as a result of aluminum complexing. Actually, seasonal 
changes in TM observations of the lakes were very small as were the 
changes in the subsurface reflectance data. The significance of these 
changes was doubtful. In addition, little seasonal change could be 
demonstrated in lake water chemistry from May to June for this data 
set. Many of these latter constituent concentrations were near the 
reported lower limit of detection. During the winter of 1986 and 
1987, the precipitation was particularly anomalous. Lack of snow 
during the winter left water levels down an average of three to four 
feet in the Sudbury area during spring, 1987. The lack of snow and 
subsequent runoff may explain the absence of a seasonal change in TM 
reflectance. More extensive seasonal observations are necessary to 
validate the season transparency hypothesis. 

An historical TM scene pair (1985-1986), however, did demonstrate 
multi-year changes that were consistent with expected changes in water 
chemistry, but lacks the chemistry and in situ optical data needed for 
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hypothesis validation. Lakes displaying the greatest TM changes are 
also the ones which were identified to be in acid sensitive strata. 

We conclude that there is likely some seasonal changes in transparency 
which can be related to the acidification process but it is also 
likely that year to year variability is significant. Strong 
relationships were found between chemical and optical properties of 
sampled lakes and the eco-physical strata within a single date. 

Optical transparency in clear acidified lakes is sensitive to water 
quality changes. 

Results show that a remote sensor can discriminate clear acid 
lakes from colored high DOC lakes based upon reflection. The clear 
acid lakes may be naturally clear. TM signals were found to be gen- 
erally higher for these lakes due to higher volume reflectance and 
greater effective transparency. Subsurface and airborne spectral 
reflectance measurements confirm this result. High DOC lakes in the 
same sensitive environments are less prone to pH change and certainly 
to changes in reflectance. Many of these lakes were originally acidic 
and will remain so but seem to be less impacted by acid deposition 
than the clearer low DOC lakes. Both lake types can be distinguished 
by remote sensing but it is necessary to first stratify the region to 
identify the acid sensitive environments. When stratification of eco- 
physical properties is used to identify acid sensitive areas TM can be 
used to pick lakes which are likely to be most sensitive to acid 
deposition and which also are indicators of temporal change. 

The opportunities for using TM to monitor multitemporal lake 
reflectance changes remains positive but additional data collections 
are considered necessary to confirm or deny the interpretations made 
in the present study. However, it is apparent that remote sensing of 
lake reflectance provides a means to identify many of these lakes and 
to possibly monitor their decline or recovery over extended period of 
time. 
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2.0 INTRODUCTION 

2.1 STATEMENT OF THE PROBLEM 

The acidification of lake waters from airborne pollutants is of 
continental proportions both in North America and Europe. A major 
problem with acid deposition is the cumulative ecosystem damage to 
lakes and forests. The number of lakes affected by this in north- 
eastern United States and on the Canadian Shield is thought to be 
enormous . 

2.2 STATEMENT OF THE OBJECTIVES 

This research had three principal objectives. First, determine 
how lake constituent concentration and lake transparency are related 
to annual acidic load. Second, investigate the utility of Thematic 
Mapper (TM) based observations to measure changes in the optical 
transparency in acid lakes. Third, examine the relationships between 
variations in lake acidification and eco-physical properties. 

2.3 BACKGROUND 

Previous investigations have suggested that DOC, which originates 
from the dissolution of humic substances, controls transparency in 
many Canadian Shield Lakes (Howard and Perley, 1982). It has also 
been established that aluminum, which is abundant in the local rocks 
and soils, is easily mobilized by acidic components contained in 
spring runoff (Hendry and Brezonik, 1984). The presence of any sig- 
nificant amount of aluminum induces a loss of DOC from the water 
column by coagulation and complexing resulting in increased optical 
transparency. This process has not been observed in lakes with normal 
pH levels associated with buffered geologies. In a normal lake, 
transparency would tend to decrease in time with the seasonal phyto- 
plankton productivity cycle. Thus seasonal changes in the optical 
transparency of lakes should potentially provide an indication of the 
stress due to acid deposition. 
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The potential for this optical response is related to a number of 
local eco-physical features with soil/geology being, perhaps, the most 
important. Other important factors include sulfate deposition, vege- 
tation type, vegetation cover, and topographic relief. The area of 
northern Ontario under study contains a wide variety of geologies from 
acid-sensitive quartzite to acid-insensitive dolomite. Annual sulfate 
deposition ranges from 1.0 to 4.0 grams per square meter (Environ- 
mental '82 Committee, 1982). 

An acidifying lake undergoes a process of decay known as oligo- 
trophication. Fewer and fewer ions of acid within the lake can be 
neutralized by the biological community. Increasing acidity further 
hampers the normal biological processes. Even though the acidity is 
not yet fatal to most fish, the lake is considered acid-sensitive and 
scientists would most like to monitor a lake at this delicate point. 

An acid-sensitive lake is thought to have, in general, high aluminum 
ion concentrations, low pH values, low alkalinity concentrations, and 
low DOC concentrations. 

Several investigators including Aimer [1974], Malley [1982], 
Schofield [1972], and Yan [1983] have reported a reduction in water 
attenuation with acidification. Aimer proposed that the changes 
resulted from probable interaction between aluminum mobilized in the 
watershed and DOC and argued that an aqueous solution with pH below 5 
will result in the precipitation of humic substances (such as DOC) 
from the water column. At pH's above 5.5 the aluminum, as aluminum 
hydroxide, will precipitate from the water column. The concentration 
of soluble aluminum will increase significantly if watershed soils are 
acidified and thus there is correlation between dissolved aluminum and 
lake pH. Acidified lakes with high concentrations of aluminum should 
also be relatively clear because of the complexing reductions of DOC. 
Aimer, however, suggests in lakes with very high humus the aluminum 
complexing does not result in precipitation. Effler's [et.al., 1985] 
description of experiments in Dart Lake not only confirm the strong 
relationship between DOC and lake transparency but also demonstrate 
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the coagulation/adsorption of DOC by aluminum. The following 
discussions relate how chemical and optical properties will be 
effected by the acidification process. 

2.3.1 PH 

Many lakes in the Northern Ontario region have experienced a 
100-fold increase in acidity (i.e., from pH=6.8 to pH=4.4) in one 
decade. Much of this is due to abnormally acidic atmospheric 
deposition and the low buffering capacity of the Shield. The present 
average acid deposition over Ontario has a pH level of 4, which is ten 
times more acidic than normal rain and 1000 times more acidic than 
neutral water. Two classifications of lakes based on pH are made most 
often. Lakes with pH's less than 6.5 are typically acid-sensitive 
lakes. These lakes have severe pH fluctuations, especially during 
spring thaw, resulting in obvious negative biotic impacts. Lakes with 
a pH of 5.0 or less can only support a few acid-insensitive plankton 
and are generally considered "acidified". Near pH 6.5 the effects are 
not as noticeable, but the pH fluctuations kill off most of the young 
biotic generations. The process leading to an "acidified" lake begins 
at a pH of 6.5. Those lakes with pH's greater than 6.5 are considered 
more or less "normal" and the water chemistry remains fairly stable 
(Environment '82 Committee, 1982). 

2.3.2 Aluminum 

Acidification transforms organic weak-acid dominated lakes to 
mineral strong-acid dominated lakes. More specifically, acidification 
decreases the availability of organic ligands for binding metals such 
as aluminum (Davis et al . , 1985). As a result, aluminum ions are 
usually found in high concentrations in acid lakes, and aluminum ion 
data could be used to predict acid-sensitive lakes. High concentra- 
tions of aluminum ions will ensure the absence of fish since aluminum 
hydroxide forms on their gills, making it difficult for the fish to 
intake oxygen. In general, if the aluminum concentrations reach 200 
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/ig/1, the lake becomes toxic to fish (Environment '82 Committee, 

1982) . 

Since precipitation has a very low aluminum concentration, the 
aluminum found in a lake's water column reflects mineral weathering 
within watersheds or mineral dissolution from lake sediments. There- 
fore, we would expect that a relationship would exist between sur- 
rounding terrain and within-lake concentrations. 

2.3.3 Dissolved Organic Carbon 

Acidified lakes found in Norway undergo a precipitation of the 
colored organic matter (DOC) in the water by acid-mobilized metals 
such as aluminum (Davis, Anderson and Berge, 1985). Increasing min- 
eral acids actually protonate organic molecules and increase their 
tendency to aggregate and precipitate. The mobilization of aluminum 
in inorganic form provides further charge neutralization of organic 
functional groups leading to their precipitation. Dissolved organic 
carbon measured from lake samples represents the amount of organics 
still within the water column and may reflect the nutrient status of 
the lake. 

2.3.4 Alkalinity 

Alkalinity is a measure of the ability of water to neutralize 
acid. The presence or absence of hydroxide, bicarbonate, and carbo- 
nate strongly influence the alkalinity or "buffering capacity" of a 
lake. Alkalinity is determined by measuring the amounts of acid 
required to neutralize alkaline water to pH 8.2 and pH 4.5 (pH 8.2 
indicates the conversion of the carbonate to bicarbonate ions and pH 
4.5 indicates the conversion of the bicarbonate ions to carbonic 
acid). These two acid levels determine the buffering capacity of the 
lake. A pH of 7.0, that of neutral water, bears little significance 
in the determination or expression of alkalinity (Chow, 1964). There- 
fore, alkalinity levels provide information not acquired with pH data 
alone. 
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When using Total Inflection Point (TIP) as a measure of alkalin- 
ity, an acidified lake is indicated when the TIP is less than or equal 
to zero (Keller and Pitblado, 1985). 

A review of the literature shows that in-lake pH levels, and con- 
centrations of DOC, aluminum and alkalinity all indicate the acid 
sensitivities of a lake. These parameters, however, are not just a 
function of in-lake processes and atmospheric loading ; they are also a 
function of terrigenous loading , i.e., a function of bedrock, soil, 
vegetation, and possibly terrain relief (Effler, Schafran, and 
Driscoll, 1985). 

2.3.5 Optical Effects 

The bio-optical state is a measure of the total effect of biolo- 
gical and chemical processes on the lake optical properties. This 
concept maintains that diverse constituents in natural waters can be 
described by a few optical parameters which represent a meaningful 
average estimate of the material present at any time and place. 

The reflectance of a lake is optically determined from the scat- 
tering and absorption processes which occur in the epilimnion (i.e. to 
the depth where the downward irradiance medium can be predicted by 
means of the radiative transfer equation). The absorption and scat- 
tering properties are inherent optical properties and do not depend on 
the light field external to the medium. There are three inherent 
properties which together are sufficient to describe the behavior of 
light in the medium. The absorption coefficient is the fraction of 
energy absorbed from the collimated beam per unit distance traversed 
in the medium. The scattering coefficient is the fraction of energy 
which is scattered out of a collimated beam per unit distance tra- 
versed by the beam. The volume scattering function describes the 
fraction of energy scattered in a specific direction per unit scatter- 
ing volume. These three inherent properties can be used to predict 
the subsurface irradiance reflectance which is described as an appar- 
ent property of the medium. The subsurface reflectance can in turn be 
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related to the above surface upwelling radiance which is also control- 
led by the radiance distribution parameters and the Fresnel transmit- 
tance. This latter radiance is a component of the radiance observed 
by an airborne radiometer or by Landsat TM. 

The scattering and absorbing agents in natural waters can be 
divided into three categories: water, dissolved materials, and sus- 

pended materials. If the absorption and scattering characteristics of 
the medium are known, the behavior of light with the suspended and 
dissolved materials in the water column can be estimated. The reflec- 
tance can be related to the constituent concentrations using a simple 
model described later in Section 7.0 since the absorption and scatter- 
ing coefficients for constituents are additive. 

For lakes in slow-weathering soil /rock conditions the amount of 
suspended mineral content is minimal. The remaining components in 
these lakes which have an optical impact are chlorophyll-a pigment and 
DOC. Both of these components have large absorption coefficients in 
the blue-green spectral region. Scattering by chlorophyll -based phy- 
toplankton is small so we are essentially dealing, in many cases, with 
an aquatic medium which is dominated by absorption. An increase in 
DOC results in increased absorption and a decrease in reflectance. 
Since the absorption cross section for DOC is large in the blue-green 
spectral region, small changes in the DOC concentration may produce 
significant changes in reflectance especially when the base concentra- 
tion is low. 

2.4 DATA COLLECTED 

Water quality parameters were measured along with in-situ optical 
data in representative lakes of the Canadian Shield. This was done to 
calibrate a Bio-Optical Model which defines the linkages between the 
acid-deposition induced chemical lake processes and the upwelling 
radiometric signals measured by the Landsat Thematic Mapper sensor. A 
spring/summer TM scene pair and companion field measurements were 
obtained for the selected study sites located in northern Ontario. 
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These data will be used to investigate possible formulations of the 
multitemporal remote sensing causal relationships between water chem- 
istry and observed changes in water transparency. 

2.5 DESCRIPTION OF THE STUDY REGION 

The study region of Northern Ontario consisted of four principal 
sites located within the following three Landsat scenes: Sudbury, 

Algoma, and Dorset. Relative locations of the study sites are shown 
in Figure 2.1 and their general characteristics are described in the 

section below. 

2.5.1 Sudbury Site 

Location : The Sudbury Site is located within the Landsat TM scene 

19-27 and has the following coordinates: 

Upper Left: 47° 40.05' -80° 49.40* 

Lower Right: 46° 16.51' -80° 36.50' 

Geology : The geology of the Sudbury site is dominated by the Lorrain 

formation which consists of quartzite, arkose, quartz sandstone, mica- 
ceous and aluminous quartz sandstone, quartz feldspar sandstone, and 
minor conglomerate and siltstone. Mafic intrusive diabase and grano- 
phyte dikes and sheets are distributed evenly throughout the site 
except near lake Wanaptei Significant amounts of conglomerate, sand- 
stone, siltstone and argillite are found in the southern half and 
northern tip of the site. In addition scattered areas of felsic intru- 
sive and metamorphic rocks, and felsic to intermediate metavolcamcs 

occur. 

Vegetation : Approximately 65% of the test site has conifer forest 

cover and approximately 35% is classified as mixed forest. 

Soil Sensitivity : Approximately 90% of this site has low potential to 

reduce acidity and the soil is predominantly shallow. The remaining 
10% of the site has a moderate potential to reduce acidity with shal- 
low soils and ultramafic bedrock. 
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Limnoloqy/Water Chemistry : The quartzite regions have very trans- 

parent lakes (e.g., Sunnywater has a Secchi depth of 25-30 meters) 
with high concentrations of aluminum, low pH values (4-5.5), low DOC 
concentrations, and metal fallout from the Sudbury smelter. The dark 
humic lakes tend to have higher pH values. 

Acid Deposition : Annual deposition in 1982 was 1.24 g/m 2 of sulfate 

2.5.2 Algoma Site 

Location ; The Algoma site is located within the TM scene 22-27 and 
has the following coordinates: 

Upper Left: 47° 21.5', -84° 25.8' 

Lower Right:47° 00.0', -84° 13.8' 

Geology : Granitic rock predominates (60%) in the Algoma site and is 

concentrated in the northeast and southwest corners. Approximately 
25% of the geology consists of acid to intermediate metavolcanics and 
15% is basic and undifferentiated metavolcanics. Several lakes are 
situated in greywacke-slate-arkose and grabbro formations. 

Vegetation : Hardwood forests predominate (Sugar Maple, Birch, Trembl- 

ing Aspen) with a few mixed stands in the lowland areas (White Birch, 
Black Spruce, and White Spruce). 

Soil Sensitivity : The northern half (approximately 55%) of the site 

has a high sensitivity to acid deposition with 0.25 to 1 meter soil 
depth with sandy texture and granite and associated alkalic bedrock. 
The southern corner (5%) is the same as the northern half of the site. 

A moderate potential to reduce acidity is found in the southern part 
of the test site (35%), which stems from a differing bedrock (ultra- 
mafic serpentine, non-calcareous silicic sediments and anorthosite) 

Limnoloqy/Water Chemistry : Lakes in this region are less transparent 

due to a higher DOC content. Levels of pH are typically between 5 and 

6 . 

2 

Acid deposition : Annual deposition of sulfate 1. 5-2.0 g/m 
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2.5.3 Dorset Site 

Location ; the Dorset site is located near the southern edge of TM 
scene 18-28. 

Geology : Acid intrusives occur throughout this area including gran- 

ite, syenite, granite gneiss, grantized sedimentary and volcanic 
rocks. 

Vegetation ; Predominantly hardwoods (Sugar Maple, Red Maple, Yellow 
Birch, Trembling Aspen) occur in this area. Hemlock and Eastern white 
pine are found in selected areas. 

Soil Sensitivity : The Dorset area is in the center of a large region 

of high deposition. West of Dorset there is less than 50% exposed 
bedrock and to the east 50 to 75% is exposed. 

Limnology/Water Chemistry : Lakes in this region are poorly buffered. 

DOC levels are higher and seech i depths are lower compared to the 
Sudbury area. 

Acid Depositio n: Annual deposition of sulfate 2.90 g/m 2 . 

2.5.4 Wawa Site 

Location : The Wawa site is located northeast of Wawa, Ontario near 

Michipicoten Bay. 

Geology : The northern third of the Wawa site consists of mafic meta- 

volcanics. Felsic metavolcanics occur in the southern tip of the site 
and are also interspersed with metasediments (conglomerate, greywacke, 
shale, arkose, and quartzite) near the middle of the site. 

Vegetation : This site contains large non-vegetated areas which have 

been impacted by the smelter fumes from Wawa. 

Soil Sensitivity : This area is primarily moderately sensitive to acid 

deposition. A small area of high sensitivity exists along the Maple 
River in the southern part of the Wawa plume. 
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LimnoloQv/Water Chemistry 

Lakes in this region are buffered , have higher pH values, high DOC 
levels, and relatively low transparency except in the immediate vicin- 
ity of the Wawa smelter plume where the lakes are acid and clear and 
highly contaminated with smelter waste. 

2 

Acid Deposition : Annual deposition in 1982 was 1.5 g/m 

2.6 SUPPORTING RESEARCH 

An historical water quality database, has been obtained from the 
Ministry of Environment for all of Ontario which contains many lakes 
within our proposed field sites. A second database is being acquired 
for approximately 300 lakes in the Sudbury area, many of which are 
located within the proposed sampling sites. The most important param- 
eters within this database are those which have impact on the optical 
transparency of the water. These parameters are chlorophyll pigments, 
suspended mineral particles, and dissolved organic carbon. Of these 
DOC is considered to have the greatest influence on optical properties 
in Northern Ontario. 

One obvious feature indicating a declining lake is low pH, but a 
low pH is not the only characteristic of an acidified lake. Chemical 
levels within a lake can also indicate its health. A study involving 
lake classification near Sudbury, Ontario used principal component 
analysis to show that chemical variability of acidified lakes is 
attributed to three main components: nutrient status, buffering sta- 

tus, and atmospheric deposition status (Pitblado et al . , 1980). 
Nutrient status of a lake could be indicated by levels of dissolved 
organic carbon, while buffering status could be indicated by the alka- 
linity of a lake. Atmospheric deposition status might be indicated by 
the annual rate of sulfate deposition within an area. 

Some historical data collected by John Fortescue at OGS, using the 
PROBAR/hel i copter over a portion of the Algoma site, were made avail- 
able to be analyzed with coincident limnological data. These data 
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were collected on August 22, 1984 and on September 6, 1985. Fortescue 
had attempted to used these data to separate clear and colored acidic 
and normal pH lakes within the site [Fortescue, 1986]. Since many of 
the same lakes were to be sampled during the August 1986 field work 
using the PROBAR radiometer, it seem reasonable to examine these data 
for potential relationships between the PROBAR measurements in TM 
bands and the measured values of DOC, pH, etc. The data set consisted 
of 113 sample locations and a representative subset was selected for 
data reduction. The reported reflectances at 10 nm intervals were 
first reduced to simulate TM band reflectances in bands 1 through 4. 
These data were then statistically correlated to the available limno- 
logical data. 

Attempts to run analyses on the combined 1984/1985 data set 
yielded very poor correlations. The 1985 data were found to be sus- 
pect because of reported instrumentation problems and further analysis 
of the 1985 PROBAR data set was therefore discontinued. The pH values 
of the 1984 data set ranged from 4.9 to 5.57 with a mean value of 
5.24. DOC values were high and ranged from 3.1 to 14.1 mg/1 with a 
mean value of 6.7 mg/1. Correlations with estimated TM reflectance 
values were considered modest (-0.73 for pH and TM band 3, -0.71 for 
pH and TM band 4). Similarly, coefficients of 0.62 and 0.64 were 
determined between the two TM bands and measured DOC. Correlations of 
comparable magnitude were observed between pH, DOC, and Seech i depth 
transparency. The lack of strong correlation was attributed to the 
relatively high levels of DOC which almost completely absorb the radi- 
ation in TM bands 1 and 2. 

2.7 STUDY ORGANIZATION 

This study was divided it into four types of activities: 1) stra- 
tification of eco-physical sensitivity, 2) water quality measurements, 

3) lake optical measurements, and 4) remote sensing measurements. 

These activities in turn supported calibration of an optical model 
which would describe the reflectance sensitivity to changes in water 
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parameters and relationships between spatial eco-physical features. 
These eco-physical features describe the environmental sensitivity to 
acidification. Our approach is outlined with the organizational flow 
chart contained in Figure 2.2. The desired result from this effort 
was to be able to identify which environments contain lakes which are 
sensitive to acidification and can be monitored using Landsat TM data 


2.8 STUDY PARTICIPANTS 

A cooperative program with Canadian agencies and Universities 
interested in the remote sensing aspects of the acid deposition prob- 
lem have resulted in an informal joint program which includes four 
major Canadian participants. These are Professor Roger Pitblado of 
Laurentian University in Sudbury, Ontario, Dr. John Fortescue of the 
Ontario Geological Survey (OGS) , Dr. Vernon Singroy of the Ontario 
Centre for Remote Sensing (OCRS), and Professor Michael Dickman from 
Brock University in Saint Catherine, Ontario. 

The Canadians are funded through the Ministry of Environment (MOE) 
and the Ontario Geological Survey for a one year period to work col- 
laboratively on the program. These funds were budgeted to support 
equally remote sensing data collection and analysis and a geochemical 
survey. 

The Canadian effort was based on meeting two separate but highly 
complementary objectives. The OGS objective was designed to look the 
relationships between environmental and geochemical studies involving 
lake acidification and remote sensing. The geochemical survey tech- 
niques developed by John Fortescue of the OGS involve analysis of 
chemical constituents in lake water samples and in bottom sediment 
cores. A mineral resource appraisal was a specific objective of the 
OGS. The MOE support was directed at examining the role remote sens- 
ing can play in the study of lake acidification in both the short and 
in the long term. The MOE had stressed that effort be placed on the 
Sudbury site where there exists an extensive limnological database. 
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The MOE plan includes examination of several historical Landsat TM and 


MSS collections. 
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3.0 ECO-PHYSICAL CHARACTERIZATION 


3.1 OBJECTIVE 

The objective of the eco-physical stratification and characteriza- 
tion of acid-sensitive parameters was to reveal the location and co- 
occurrence of environmental attributes that influence lake acidifica- 
tion. The study areas were stratified into the following four 
parameters : 

1. type and percent cover of vegetation, 

2. soil and bedrock buffering capacity, 

3. topographic relief, 

4. sulfate deposition rate. 

The acid sensitivities of these areas were then determined, based 
on these four parameters. Each of these parameters affects the sensi- 
tivity of the ecosystem a lake is found in and ultimately affects the 
water chemistry and optical signature of that lake. Stratification 
also provided a basis to characterize lakes within study areas which 
aided in the sampling design. 

3.2 PROCEDURE 

The three Landsat scenes were stratified into eco-physical units, 
or "polygons", based upon soil/bedrock sensitivity, vegetation sensi- 
tivity, topographic-relief sensitivity and acid- deposition sensitiv- 
ity. Sensitivity values were assigned to each polygon and combined in 
a linear function which produced a "sensitivity index" for each poly- 
gon using a sensitivity model. Max imum- 1 ikel i hood clustering of these 
sensitivity indexes then revealed the location and co-occurrence of 
similar polygons. 

3.3 STRATIFICATION OF ECO-PHYSICAL FEATURES 

The Algoma, Sudbury, and Dorset study areas were stratified in 
terms of bedrock/soil, vegetation, relief and sulfate deposition. 
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Four mylar overlays were constructed, one for each of the variables, 
at a scale of 1:250,000. 

3.3.1 Vegetation and Percent Cover 

The lowest pH values are found in coniferous forests. Fir trees 
are often found growing on weathering-resistant soils and bedrock. 

When precipitation falls on this type of area, the acidic water flows 
largely unaltered into nearby lakes at a pH of 5.6. Broadleaf forests 
are generally found in terrain of higher pH, so precipitation is neu- 
tralized more before it enters a lake. A much higher rate of sulfate 
deposition would be necessary to make the pH of runoff from a decid- 
uous forest reach that of a coniferous forest (Environment '82 Com- 
mittee, 1982). 

Percent cover of vegetation also plays a factor in lake acidifica- 
tion. If percent cover is low, the extent and volume of surface run- 
off is frequently higher than for average cover conditions increases. 
Under these conditions, very little of the precipitation has time to 
penetrate into the rock and/or soil and become neutralized by the 
buffering systems. 

TM satellite images were used for vegetation classification and 
lines were drawn between areas of different vegetation types and dif- 
ferent percent covers of these types. Vegetation was categorized as 
conifer, hardwood, mixed or barren. If an area's vegetation consisted 
of 80% or more of either conifer forest or hardwood forest, then it 
was classified hardwood or conifer, otherwise it was classified as a 
mixed forest. 

Percent cover for an area was derived using existing soil and 
bedrock sensitivity maps published by the Environment Canada Lands 
Directorate in 1983. These maps outline percent exposed bedrock at 
three levels: 0-24%, 25-50%, and 50-99%. Since there were no exten- 

sive areas of low vegetation, such as prairies, marshes, etc., the 
following equation was used: 

(Percent forest cover) = 1 - (Percent exposed bedrock) . 
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Percent forest cover was divided into three classifications: 

1. 0 - 49 % cover, 

2. 50 - 74 % cover, 

3. 75 - 99 % cover. 

Vegetation and percent cover sensitivities were derived from the lit 
erature (Environment '82 Committee) and are shown in Table 3.1. 

TABLE 3.1. VEGETATION AND PERCENT COVER SENSITIVITIES 

Cover 

hardwood 
hardwood 
hardwood 

mixed 
mixed 
mi xed 

conifer 
conifer 
conifer 

These sensitivity values rank the combinations of vegetation type and 
percent cover on a scale from 1 to 10. Terrain with conifer forest 
cover was rated most sensitive and terrain with hardwood forest cover 
was rated least sensitive. The higher the percent cover the less 
sensitive the polygon was rated for potential damage. 

3.3.2 Sulfate Deposition 

Large emissions of sulfur dioxide and nitrogen oxide from combus- 
tion (usually within coal burning industries) lead to their oxidation 
in the atmosphere to sulfuric acid and nitric acid. These acids dis- 
solve in water droplets and fall to the ground via some form of pre- 
cipitation. The presence of sulfuric acid in precipitation over the 
Continental Shield results in 100 times more acid entering these 
already poorly buffered ecosystems (Hendry and Brezonick, 1984). 
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The sulfate deposition overlay was drawn from enlarged 1981 
meteorologic maps (Chan, et al. 1983) provided by the Ontario Ministry 
of the Environment (see Figure 3.1) Sulfate deposition was measured 
in grams/m 2 /year. Across all three areas, the following six classifi- 
cations were derived from the maps in terms of deposition rates: 

1 . 1 . 0 - 1 . 5 , 

2. 1. 5-2.0, 

3. 2. 0-2. 5, 

4. 2. 5-3.0, 

5. 3. 0-3. 5, 

6. 3. 5-4.0. 

Sulfate deposition was assigned sensitivity values based on amount 
of sulfate deposited. Each of the six levels was assigned equally 
spaced sensitivity values on a scale from 1 to 10. The highest sul- 
fate deposition was given the highest sensitivity value. The results 
are given below in Table 3.2. 


TABLE 3.2. 


SENSITIVITY VALUES OF SULFATE DEPOSITION LEVELS 


gm/nr/year 

Sensitivity 

1.0-1. 5 

1.67 

1. 5-2.0 

3.33 

2.0-2. 5 

5.00 

2. 5-3.0 

6.67 

3. 0-3. 5 

8.33 

3. 5-4.0 

10.00 


3.3.3 Bedrock and Soil 

In general, the easier the ground materials around a lake weather, 
the less susceptible that lake is to acidification. Thus, weather- 
ability of the lake's surrounding bedrock and soil play a large factor 
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Figure 3.1 . The Annual Deposition (G/M**2) of Sulfate in Ontario 
(from Chan, Tang and Lusis, 1983). 
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on the lake's acidity. The rate at which bedrock and soil weather 
depend on their hardness and their ability to release buffering ions 
which counter lake acidification by reducing the impact of the water 
runoff. 

Bedrock resistant to weathering does not neutralize acid rainwater 
therefore it is associated with acidic lake systems. Sensitivities 
for bedrock/soil combinations were derived from the Environment Canada 
Sensitivity Maps. Bedrock was divided into four categories based on 
its sensitivity. These four categories are found in Table 3.3. 


M 


TABLE 3.3 BEDROCK SENSITIVITY CATEGORIES 
Description 

limestone, marble, dolomite 

carbonate-rich siliceous sedimentary: shale, limestone; 

noncalcareous siliceous with carbonate interbeds: shale, 

siltstone, dolomite; quartzose sandstone with carbonates. 

ultramafic rocks, serpentine, noncalcareous siliceous 
sedimentary rocks: black shale, slate, chert; gabbro, 

anorthosite: gabbro, diorite; basaltic and associated 

sedimentary: mafic volcanic rocks. 

granite, gneiss, quartzose sandstone, syenitic and asso- 
ciated alkalic rocks. 


The ability of the soil to neutralize the acid was found to be 
the most important factor influencing the susceptibility of a lake to 
acidification. Lime-rich, easy-weathering soils protected the lakes, 
but lakes surrounded with sandy soil and expanses of flat bare rock 
are mostly acid (Environment '82 Committee, 1982). Basically three 
categories of soil can be defined: easy-weathering clay, normal - 

weathering loam, and resistant-weathering sand. 

The soil's depth also affects the neutralization of precipitation. 
A deeper soil will contain larger quantities of weatherable minerals 
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and other buffering substances. Thin soils are often leached of such 
buffering substances. In the stratification, one of the soil types 
(clay, loam or sand) was assigned to each polygon. Each polygon was 
also assigned a unique soil depth. The soil depth categories used are 
shown in Table 3.4. 

TABLE 3.4. SOIL DEPTH CATEGORIES 
Definition 

> 1 m average soil thickness 
25 cm - 1 m average soil thickness 
< 25 cm average soil thickness 

Different combinations of bedrock type, soil type, and soil depth were 
already ranked on the Environment Canada maps from most to least sen- 
sitive. Since there were 28 soil /bedrock combinations, the most sen- 
sitive combination was assigned a 10.0. The other combinations were 
assigned sensitivities ranging from 1 to 10 separated by units of 
10/28. These combinations are shown in Table 3.5. 
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TABLE 3.5. BEDROCK/SOIL SENSITIVITY INDEX VALUES 


ROCK TYPE 

SOIL TYPE 

SOIL DEPTH 

SENSITIVITY VALUE 

1 

clay 

deep 

.36 

1 

loam 

deep 

.71 

1 

sand 

deep 

1.07 

1 

clay 

shallow 

1.43 

1 

loam 

shallow 

1.79 

1 

sand 

shallow 

2.14 

1 

none 

bare 

2.5 

2 

clay 

shallow 

2.86 

3 

clay 

shallow 

3.21 

2 

clay 

deep 

3.57 

3 

clay 

deep 

3.93 

4 

clay 

deep 

4.29 

2 

loam 

deep 

4.64 

3 

loam 

deep 

5. 

2 

sand 

deep 

5.36 

3 

sand 

deep 

5.71 

2 

loam 

shallow 

6.07 

3 

loam 

shal low 

6.43 

2 

sand 

shallow 

6.79 

3 

sand 

shallow 

7.14 

2 

none 

bare 

7.5 

3 

none 

bare 

7.86 

4 

clay 

shallow 

8.21 

4 

loam 

shallow 

8.57 

4 

loam 

deep 

8.93 

4 

sand 

deep 

9.29 

4 

none 

bare 

9.64 

4 

sand 

shallow 

10.00 

3.3.4 Relief 

Since the 

extent and volume of surface runoff plays an important 

factor in lake 

acidification 

, the topographic relief of the terrain 

surrounding a 

lake would help determine its 

acidification state. An 

area with steep topographic 

relief would allow less time for precipi- 

tation to penetrate the soil 

and bedrock and 

become neutralized. Flat 

topographic relief would contribute more to 

the neutralization of 

precipitation 

since the extent and volume of 

surface runoff would be 

less. 

Relief was 

divided into 

three categories 

: steep, rolling, and 

level. This information was 

extracted from 

standard topographic maps 
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at a scale of 1:250,000. Change in elevation across unit distances 
was measured perpendicular to elevation contours and categorized into 
one of three types for each polygon. These categories are shown in 
Table 3.6. 


TABLE 3.6. TOPOGRAPHIC RELIEF CATEGORIES 

Category Definition 

level: < 400 ft change in 2 kilometers 

rolling: > 400 ft < 800 ft change in 2 km 

steep: > 800 ft change in 2 kilometers 


Topographic relief levels were assigned three sensitivity values, 


Table 3.7, 


from 1 to 10 

. These three values are 

TABLE 3.7. 

RELIEF SENSITIVITY VALUES 

Relief 

Sensitivity Value 

level 

3.33 

rol ling 

6.7 

steep 

10.00 


3.4 COMPOSITE MAP CONSTRUCTION 

The four maps were produced for each of the ecosystem parameters 
(bedrock and soil, sulfate deposition, terrain relief, and vegetation 
type and percent cover). Each map consisted of polygons that repre- 
sented uniform ecosystem parameters and that were assigned correspond- 
ing sensitivity values. A composite map was then produced for each of 
the study areas by overlaying the four ecosystem parameter maps, and 
tracing them on to one overlay (see Figure 3.2). Ultimately, the new 
polygons created with the composite map had four sensitivity values: 
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Figure 3.2. The Stratification Procedure. 
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one for bedrock/soil, one for vegetation, one for relief and one for 
the sulfate deposition. 

The three composite maps produced 694 polygons with a minimum 
polygon size of 25 square kilometers. Each polygon was numbered from 
1 to 694. A computer program was written and used to read the polygon 
number, forest type, percent cover, bedrock type, soil type, soil 
depth, topographic relief and sulfate deposition into computer memory. 
A program subroutine was used to assign four ecosystem sensitivity 
values, ranging from 1 to 10, to each polygon and compute the sensi- 
tivity index for each polygon using the sensitivity index model. 

A list of the polygons with eco-physical characteristics and sen- 
sitivity index values is found in Appendix A. 

3.5 SENSITIVITY INDEX MODEL 

A sensitivity index model was developed which assigned a sensitiv- 
ity index to each composite map polygon. The sensitivity index, SI, 
is a function of a linear combination of the four ecosystem parameters 
within the polygon: 

SI = A x (bedrock/soil sensitivity value) 

B x (vegetation sensitivity value) 

C x (sulfate deposition sensitivity value) 

D x (topographic relief sensitivity value). 

The coefficients A, B, C and D were derived from the literature, but 
in the absence of quantitative information. An ecosystem sensitivity 
study in Sweden concluded that bedrock and soil were found to be the 
most important factors influencing the susceptibility of a lake to 
acidification (Environment '82 Committee). Also, areas of nearly 
equal rates of sulfate deposition, but differing types of bedrock and 
soil, have been found to contain lakes of different buffering capaci- 
ties, supporting the idea that bedrock and soil are the most important 
eco-physical parameters in terms of lake sensitivity. Therefore the 
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coefficient "A" equals four, the highest number assigned to a coeffi- 
cient. A review of the literature indicated that vegetation type was 
highly correlated with soil and bedrock type in terms of sensitivity, 
so the vegetation sensitivity value was weighted as the second most 
important variable. 

If the vegetation and soil/bedrock sensitivity values were iden- 
tical in two areas, it is assumed that sulfate deposition would affect 
the sensitivity of a lake within the area more than topographic relief 
would. Therefore the following equation was developed: 

SI = 4 x (bedrock/soil sensitivity value) 

3 x (vegetation sensitivity value) 

2 x (sulfate deposition sensitivity value) 

1 x (topographic relief sensitivity value). 

The sensitivity index of an eco-physical polygon is driven by the 
bedrock/soil and vegetation sensitivity values. The sulfate deposi- 
tion and topographic relief sensitivity values still contribute to an 
area's sensitivity, so they are included in the model but weighted as 
less important. Therefore, it is hypothesized that the sensitivity 
index rates the acid sensitivity of an eco-physical area on a scale 
from 1 to 10. 

3.6 CLUSTERING OF MODEL SENSITIVITY VALUES 

The sensitivity indexes of the polygons (approximately 694) were 
then clustered using a maximum likelihood hierarchical clustering 
procedure. The results of this clustering procedure has produced 10 
significantly (p > .95) different clusters (see Appendix B) . These 
clusters are summarized in Table 3.8. 
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TABLE 3.8 SENSITIVITY RATINGS AND TYPE VALUES FOR THE TEN 
SIGNIFICANTLY DIFFERENT CLUSTERS 


CLUSTER RATING 

BEDROCK/SOIL 

VEGETATION 

RELIEF 

SULFATE 

DEPOSITION 

1 

5.66 

7.04 

4.67 

5.57 

4.40 

2 

6.36 

8.05 

4.65 

5.78 

5.82 

3 

6.74 

8.16 

5.83 

5.28 

5.00 

4 

6.02 

7.67 

4.63 

5.25 

5.18 

5 

7.41 

8.47 

7.13 

5.62 

6.59 

6 

3.55 

3.28 

2.08 

5.57 

5.27 

7 

7.07 

8.50 

6.37 

5.36 

6.10 

8 

5.14 

5.96 

4.71 

5.46 

3.97 

9 

7.83 

8.71 

8.53 

5.20 

6.29 

10 

4.34 

5.21 

3.82 

5.01 

3.05 


The ten clusters are described in terms of their mean eco-physical 
sensitivity values in the following paragraphs. 

Cluster 1 is characterized by shallow sandy soils over rock types 3 
and 4 with less than 50% cropping out. Vegetative cover is a mixture 
of conifers and hardwoods with a dominance of the hardwoods. The 
terrain is level to rolling. The average acid deposition is approxi- 
mately 2.0 g/m /yr. 

Cluster 2 is characterized by moderate depth soils over rock type 4 
with less than 50% cropping out. Vegetative cover is a mixture of 
conifers and hardwoods with a dominance of the hardwoods. The terrain 
is level to rolling. The average acid deposition is approximately 2.5 
g/m 2 /yr. 

Cluster 3 is characterized by deep sandy soils over rock type 4 with 
less than 50% cropping out. Vegetative cover is a mixture of conifers 
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and hardwoods with a dominance of the hardwoods. The terrain is level 

p 

to rolling. The average acid deposition is approximately 2.5 g/m /yr. 

Cluster 4 is characterized by moderately deep soils over rock type 4 
with less than 50% cropping out. Vegetative cover is a mixture of 
conifers and hardwoods with a dominance of the hardwoods. The terrain 
is level to rolling. The average acid deposition is approximately 
2.25 g/m 2 /yr. 

Cluster 5 is characterized by moderately deep sandy soils over rock 
type 4 with less than 50% cropping out. Vegetative cover is a mixture 
of conifers and hardwoods with a dominance of the hardwoods. The 
terrain is level to rolling. The average acid deposition is approxi- 

O 

mately 2.75 g/nr/yr. 

Cluster 6 is characterized by deep clay soils over rock type 3 with 
less than 30% cropping out. Vegetative cover is mostly hardwood. The 
terrain is level to rolling. The average acid deposition is approxi- 

p 

mately 2.25 g/m /yr. 

Cluster 7 is characterized by shallow sandy soils over rock type 4 

with less than 50% cropping out. Vegetative cover is a mixture of 

conifers and hardwoods with a dominance of the conifers. The terrain 
is level to rolling. The average acid deposition is approximately 2.5 
g/m 2 /yr. 

Cluster 8 is characterized by moderately deep sandy soils over rock 
type 3 with less than 50% cropping out. Vegetative cover is a mixture 
of conifers and hardwoods. The terrain is level to rolling. The 

p 

average acid deposition is approximately 2.0 g/m /yr. 

Cluster 9 is characterized by shallow sandy soils over rock type 4 

with less than 25% cropping out. Vegetative cover is dominated by 

conifers. The terrain is level to rolling. The average acid deposi- 

p 

tion is approximately 2.5 g/m /yr. 

Cluster 10 is characterized by deep sandy soils over rock types 3 and 
4 with less than 50% cropping out. Vegetative cover is a mixture of 
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conifers and hardwoods with a dominance of the hardwoods. The terrain 
is level to rolling. The average acid deposition is approximately 1.5 
g/m 2 /yr. 

These clusters are separated by only small changes in the mean 
value for each sensitivity index. The standard deviations of the 
above mean sensitivity index values was typically only one or two 
percent. Each cluster was color coded as shown in Figure 3.3. Color 
coded maps that show the location of the polygons within each cluster 
are shown in Figures 3.4 3.5 and 3.6. The listing of all eco-physical 
polygons by cluster with the strata descriptors is given as Appendix 
A. The summary statistics for the clusters is given in Appendix C. 

The above clusters were further grouped into three classes which 
are shown in Table 3.9. 

TABLE 3.9. CLUSTER CLASSES 
Class Clusters 

insensitive 1, 6, 8, 10 

mildly sensitive 2, 3, 4 

sensitive 5, 7, 9 

3.7 SAMPLE SITE SELECTION 

Site selection for in situ lake measurements was based upon the 
stratification and clustering analysis described above and each of the 
following considerations: (1) availability of historical water quality 
and remote sensing data, (2) existing Canadian initiatives to collect 
site-specific data, (3) accessibility, and (4) coverage of eco- 
physical lake types. Sites selected included (1) Algoma, (2) Sudbury, 
(3) Wawa, and (4) Dorset. Nine of the ten clusters were represented 
by the selected sites. 

The Canadian program recommended the use of the Algoma and Sudbury 
sites, each comprising approximately 1000 sq. km. Priorities were set 
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Figure 3.4 The Algoma Area Clusters and Sampling Sites 


ORIGINAL PAGE 
COLOR PHOTOGRAPH 


ORIGINAL PMSe & 
CsT ^OR QU%ITY 


P»ec£Dl&Q PAGE BLANK NOT Filtm/ 




39 




OSJG'faAL s-.'A-JL. 

,, PHOTOGRAPH 



il.lklTI 


Figure 3.5 The Sudbury Area Clusters and Sampling Site 
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for each of the four collection sites based upon group interests and 
availability of resources. First priority was given to the Sudbury 
site, second to Algoma, and third to Wawa. The Dorset site was viewed 
to be largely beyond the reach of a one-month field program and would 
only be addressed after the other data objectives had all been met. A 
lake sampling budget of approximately 300 samples was divided between 
the first three sites with 150 samples allocated to Sudbury, 130 allo- 
cated to Algoma, and 20 to Wawa. An additional 25 samples would be 
taken to support the Dorset sampling if resources were available. 


RECEDING PAGE BLANK NOT FILMED 
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4.1 LAKE SAMPLING STRATEGY 

The ERIM field plan specified sampling at three different levels 
and with three different optical measurements. Field data collections 
were made during the summer of 1986 and spring of 1987. The August 
1986 collections included three sites: Sudbury, Algoma, and Wawa. At 
each of these sites, water samples were gathered from a well distrib- 
uted set of lakes using a helicopter. Radiometric measurements were 
made using Landsat TM, a helicopter (BELL-206) spectral radiometer 
(PROBAR), a subsurface spectral irradiance meter, and a subsurface 
beam transmissometer. The sampling strategy was to gather subsurface 
measurements from a small number of lakes and in sufficient number to 
calibrate a subsurface reflectance model. Airborne spectral measure- 
ments were gathered over a much larger set to be used to extend the 
subsurface results to a broader set of lake conditions. Finally these 
lake reflectance spectral characteristics were used to predict the 
reflectance characteristics of the still larger TM lake sample data 
set. The strategy in this three-tier sampling scheme was to develop a 
model/relationship from the in situ optical measurements and the mea- 
sured limnological parameters. This "optical response model " , once 
validated, was extended to the PROBAR data set and finally to the 
Landsat data set where it aided in the interpretation of TM 
observations. 

During August 1986 field data were gathered from each of the three 
sites which included 21 water quality parameters (296 lakes), detailed 
subsurface optical measurements (12 lakes), airborne spectral radio- 
meter measurements (102 lakes), and Landsat data. Most of these mea- 
surements were made in the Algoma and Sudbury sites (shown as Figures 
D.l and D.2). All water chemistry data are compiled as Appendix D. 
PROBAR spectral radiometer measurements were made in most of the lakes 
that were larger than 20 hectares. The subsurface optical measure- 
ments were made in a representative set of lakes at each site. Water 
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parameters were determined from collected samples by the MOE on-site 
or at the Toronto Laboratory. Water parameters especially important 
to this study included DOC, conductivity, total chlorophyll-a pigment 
concentration, pH, sulfate, alkalinity, TIP, turbidity, suspended 
solids, and aluminum. 

The May-June 1987 field effort involved collecting subsurface MER 
reflectance and transmissometer data on four separate dates from eight 
lakes. Water samples were also collected and were processed by the 
MOE. Field data collections were made on 5 May, 14 May, 13 June, and 
29 June at four to eight lakes in the Sudbury site. These data were 
collected coincident with the TM overpass on each of those dates. Two 
of these TM acquisitions (12 May and 13 June) were of excellent qual- 
ity and were requested from NASA GSFC. No PROBAR airborne radiometer 
data were collected during the spring period because the unit was not 
available for project use. 

4.2 SUBSURFACE OPTICAL MEASUREMENTS 

Two instruments were used to make the subsurface optical measure- 
ments: a subsurface spectroradiometer (Biospherical Inc. MER-1000) 
with 11 narrow spectral bands (410, 441, 488, 520, 540, 560, 589, 625, 
671 and 694 nm) and a transmissometer (SEATECH Inc.) with a single 
wavelength at 664 nm. These instruments were used to characterize the 
optical properties in several of the PROBAR-sampled lakes. 

The MER-1000 subsurface upwelling and downwelling spectral spec- 
tral scans were collected in the field at variable sampling depths 
below the lake-water surface. MER data collections were made from a 
canoe (August 1986) and from a float plane pontoon (May-June 1987). 

The canoe measurements each consisted of 20 scans and the float plane 
measurements consisted of 10 scans. Fewer scans were used during the 
plane measurements since the instrument was allowed to drop through 
the water column at a faster rate. At each station a series of 
upwelling and downwelling irradiance measurements were made in suc- 
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cession. A pressure sensor in the MER recorded the depth of each 
spectral scan. 

4.3 AIRBORNE RADIOMETER MEASUREMENTS 

A helicopter-mounted (BELL 206) spectroradiometer (PROBAR) was 
used to collect radiometric data in each of 10 narrow spectral bands 
(443, 470, 520, 550, 580, 610, 640, 670, 700 and 732 nm) at the center 
of each sample lake. 

PROBAR data was collected on four days in 1986: 

August 12 15 Lakes 

August 13 54 Lakes 

August 14 18 Lakes 

August 18 46 Lakes 

Lakes sampled with the PR0BAR were limited to those large enough 
to be visible in TM imagery and sufficiently deep not to produce a 
bottom reflected signal. The PR0BAR unit had been rented from Moniteq 
Ltd., Toronto, Ontario and was controlled with an IBM PC that also 
was mounted in the helicopter. The PC logged the radiometer data and 
allowed easy transfer to the DEC VAX780 for data analysis. 

4.4 LANDSAT TM ACQUISITIONS 

All possible Landsat TM acquisitions were requested for the 
Algoma, Sudbury, and Dorset scenes for the month of August 1986. 

Algoma and Sudbury coverage were requested for May and June 1987. Of 
the scenes collected, four were considered sufficiently cloud- free to 
be useful. Image tapes were obtained from NASA GSFC Landsat office 
and are listed in Table 4.1. 
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TABLE 4.1. IMAGE TAPES REQUESTED FROM NASA GSFC LANDSAT OFFICE 


Path/Row 

Date 

19/27 

August 13, 1986 

19/27 

May 12, 1987 

19/27 

June 13, 1987 

22/28 

August 18, 1986 


All of the other acquisitions were considered non-usable based upon 
the positive print of TM band one received from GSFC. 

4.5 DATA QUALITY MEASURES 

Provisions were made to ensure the quality of the data measure- 
ments. During the MER data collection, deck cell measurements of 
downwelling hemispherical irradiance were taken coincidentally. This 
ensured that the MER downwelling and upwelling profile measurements 
were taken while the downwelling irradiance remained constant. 

When TM signals were being extracted, band four signals of water 
surfaces were examined for high standard deviations (> 0.5). If the 
standard deviation was higher than 0.5, it was assumed that the data 
were contaminated with either bottom or land reflectance, and they 
were not used. 

Before transmissometer measurements were made, the air voltage was 
checked and recorded. The transmissometer measurement was only made 
if the air voltage was in the appropriate range. This air voltage was 
later used for calibration when calculating attenuation coefficients. 

PROBAR measurements were corrected for the time of day and were 
calibrated using a white card of known reflectance. Instrument cali- 
bration was also done in the lab before the field work. 
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5.0 SUBSURFACE AND AIRBORNE RADIOMETRIC DATA REDUCTION 

Radiometric data collected with the Biospherical MER-1000 radiom- 
eter, the SeaTECH transmissometer, the PROBAR spectral radiometer, and 
Landsat TM were reduced as described in the following sections. 

5.1 MER DATA REDUCTION 

MER-1000 data were first used to interpolate the irradiance data 
to common depths on a logarithmic scale before computing values of 
subsurface reflectance. The slope of the depth log-irradiance regres- 
sion equation defines the average irradiance attenuation coefficient 
(K). The irradiance attenuation coefficient changes very little 
within the mixed layer, but rapidly within the transition zone (ther- 
mocline). The thickness of the mixed layer was easily determined from 
the temperature depth profile. Therefore only irradiance measurements 
from the mixed layer were used to determine K. Downwelling irradiance 
attenuation for low DOC lakes (Sunnywater and Wolf) and high DOC lakes 
(Whitepine and Barbara) are shown in Figure 5.1. 

Subsurface spectral reflectances were calculated at 2, 4, 6, and 8 
meters below the surface. Example reflectance curves are shown in 
Figure 5.2, along with the DOC and Chlorophyll-a measurements. The 
impact of DOC and Chlorophyll-a on reflectance is apparent. As DOC 
increases the blue-green portion of the reflectance spectrum is dimin- 
ished due to highly selective absorption. Chlorophyll-a also dimin- 
ishes the measured reflectance below 520 nm, due to absorption. Wave- 
lengths greater than 520 nm absorption are reduced and backscattering 
is increased. The reflectance calculations at 700 nm are not con- 
sidered valid since the irradiances are very small and contaminated by 
sensor noise. 

In the spring of 1987 the MER pressure sensor was calibrated so 
measurement depths were available without depth correction. The pres- 
sure sensor in August 1986 sampling period was precise but it was not 
accurate. A control profile was made during which actual and measured 
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Figure 5.2 Subsurface Reflectance R (X) 
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depths were recorded and a simple linear relationship was found 
between them. 

Depth = Measured Depth / .567 

To obtain reflectance values it was necessary to develop two 
linear equations describing the relationship between the natural 
logarithm of irradiance (In (E) ) and corrected depth for both the 
upwelling and downwelling profiles. The diffuse attenuation coeffi- 
cient determines the rate of irradiance loss through the water column 
and is defined by the following equation. 

K <» -tihyf 

The irradiance data collect at multiple depths were first used to 
estimate K from the solution to the above equation as given by the 
following linear form. 

ln(E(X,z)) = K(X)*z + intercept 

Depths of 2, 4, 6 and 8 meters were then entered into the linear 
equation to estimate In (E u ) and ln(E d ). Reflectance at these four 
depths were then produced using the following equation: 

R(X,z) = EXP(ln(E u (X f z)) - ln(E d (X,z)) 

Where E u (X,z) = upwelling irradiance at z meters 
and Ed(X,z) = downwelling irradiance at z meters. 

5.2 TRANSMISSOMETER DATA REDUCTION 

SeaTECH transmissometer profiles were made at every station coin- 
cident with the MER measurements. Voltage measurements were made 
usually at 2, 4, 6, and 8 meters after an air reading was made at each 
station. 

Corrected voltage was then obtained using the following equation: 
Cvolt = (Lab Air/Field Air) x (Mvolt - .003) 
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where Cvolt = Corrected voltage 

Lab Air = Lab air reading = 4.775 volts 
Field Air * Field air reading 
Mvolt = Measured voltage 

Fractional transmission could be determined since it is known that 
100% transmission through 25 cm of pure water has a corrected voltage 
of 5 volts. Fractional transmission through 25 cm of lake water is 
found using the following relationship: 

T(664nm) = (Cvolt) / 5 volts. 

The beam attenuation coefficient (c) can be derived using the frac- 
tional transmission in the following equation: 

c(664nm) = -4 L0G(T(664nm) 

The reduced transmission and beam attenuation coefficients for all 
SeaTECH measurements are given as Appendix E. 

5.3 PROBAR DATA REDUCTION 

One objective in reducing the PROBAR data was to estimate illumi- 
nation independent reflectance values which could be compared to the 
MER data derived values. The airborne PROBAR measurements, however, 
were made complicated by the helicopter blade motion and by the need 
for irradiance reflectance given the PROBAR is a radiance device. The 
rotating blade interfered with the downwelling irradiance meter and 
also possibly with the upwelling radiance measurements as well. The 
raw data from several dates showed a significant change in downwelling 
irradiance between measurements taken on the ground using a standard- 
ized white reflectance card. This effect was dependent on time of day 
and date illumination conditions. These conditions necessitated a 
series of five corrections be made to these data in order to make them 
compatible to the MER reflectance data. These corrections were (1) 
for standardized white card reflectance, (2) for airborne conditions, 
(3) for time of day, (4) for day-to-day variations in sky illumina- 
tion, and (5) for surface reflectance. 
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Upwelling radiance, L U (X), and downwelling irradiance values were 
read for ten 20 nm - wide bands ranging from 433 nm to 710 nm. 
Reflectance was computed in the following manner: 

R(X,0) = M(X,0)/E d (X,0) 
where M(X,0)= L u (X)*ir 

All dates show a large change in downwelling irradiance between 
measurements taken on the ground (white card measurements) and mea- 
surements taken when the helicopter was airborne (all lake measure- 
ments). This discrepancy was accounted for in the change in heli- 
copter blade tilt. When the instrument was airborne, the blades were 
tilted at a higher angle, thus allowing more light to reach the down- 
welling irradiance sensor. A correction was made by producing a 
second order regression equation of all airborne downwelling irradi- 
ances as a function of time. The true white card downwelling irra- 
diance was then estimated using the resultant equation. This correc- 
tion was made for each PROBAR band. 

All data needed to be normalized to one unique white card reflec- 
tance for each band. The white card used for correcting the data was 
known to have a nearly constant reflectance value (.989) for the bands 
being studied. The white card reflectances were fit to a second order 
equation using time as the independent variable producing the measured 
white card reflectance curve. The true lake reflectance is adjusted 
by the same percent difference as that between the measured white card 
reflectance (MWCR) the known white card reflectance curve. 

R(true) = R (measured) x [l - [ ^^"r* 989 ]] 

A final correction was made to PROBAR measurements which was lake- 
dependent. The assumption was made that no internal lake reflectance 
was measured in the band centered at 700 nm. This measurement was 
assumed to be an indication of wave induced surface reflected noise 
and thus was subtracted at all wavelengths. This correction only 
changed the offset of the spectral reflectance curve, not its shape. 
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The above and below surface corrected PROBAR reflectances are given as 
Appendix B. 
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6.0 LANDSAT TM PROCESSING METHODS 
6.1 LAKE SIGNATURE EXTRACTION 

Extraction software was applied to all three TM scenes. Lake 
signals were extracted from the TM images by finding the latitude and 
longitude of lakes of interest on topographic maps and using these 
latitudes and longitude to extract lake signatures from geometrically 
corrected imagery using extraction software. Nine brightness values 
were extracted from each lake and their means were used in subsequent 
processing. A three by three pixel area was extracted and the mean 
signal and its standard deviation for each band were recorded. To 
ensure that the spectral signatures represented water and not cloud, 
shoreline or bottom reflectance, TM band 4 signals were inspected. 
Average signals in TM band four were found to range between 11.0 and 
14.0 with a standard deviation for values within an individual samples 
of less than 1.0. Thus for samples which had mean values outside this 
range or with sample standard deviations greater than 1.0 the sample 
was rejected and considered to indicate a non-water mixed reflectance. 
The rejected samples were replaced with values extracted from another 
part of the lake surface. Brightness values were extracted from the 
approximate center of each lake based upon the latitude and longitude 
of each lake center. These extracted mean values were then correlated 
to historical water chemistry data available for the same lakes as 
discussed in Section 8.0. 

The TM data extracted is summarized in Table 6.1. 


TABLE 6.1. 

THEMATIC MAPPER 

DATA EXTRACTED 

Path/Row 

Quad 

Date 

22/27 

1 

8/18/86 

22/27 

4 

8/18/86 

22/27 

4 

5/27/85 

19/27 

3 

8/13/86 

19/27 

3 

5/22/85 
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6.2 SOLAR ELEVATION ANGLE CORRECTION 

All lake data were corrected for the solar elevation angle of each 
scene. This correction simply involved dividing each brightness value 
mean by the cosine of the solar zenith angle. 

6.3 ATMOSPHERIC HAZE CORRECTIONS 

A haze correction needed to be applied to the TM data so that real 
comparisons could be made between lakes within and between scenes 
which had varying amounts of haze distorting the signals. Lakes of 
equivalent Dissolved Organic Carbon (DOC) concentrations should have 
similar TM signals in band one but these data showed instead wide 
variations. The lakes with elevated TM band one counts also had ele- 
vated counts in bands two, three, and four. Since band 4 counts 
represent virtually no internal lake reflectance, it was hypothesized 
that relative differences between lakes in band four represented dif- 
ferences in atmospheric haze. Linear regression analyses between 
bands one and four, bands two and four, and bands three and four 
showed nearly linear behavior but with different slope and a small 
intercept. Also, these derived slope values were found to be scene 
dependent. The slopes between bands were derived using regression 
analyses and used directly in the haze correction algorithms. Thus 
the correction for haze was both wavelength dependent and scene depen- 
dent. The following three equations are the haze correction algo- 
rithms for the three TM bands used: 

TM-1 (com) = TM-1 - ( TM-4 x Mj) 

TM-2(corr) = TM-2 - ( TM-4 x M 2 ) 

TM 3 (corr) = TM 3 - ( TM 4 x M3) 

M l* M 2 , and M3 are the slopes between bands one and four, bands two 
and four, and bands three and four, respectively. 

This procedure reduced the impact of haze as indicated by the 
improved correlation between TM band one signals and DOC (i.e. from 
0.62 to 0.83). 
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7.0 DEVELOPMENT OF A BIO-OPTICAL REFLECTANCE MODEL 
7.1 REFLECTANCE MODEL 

A TM radiative transfer model was developed to predict possible 
changes in radiometer signal levels which result from field-measured 
changes in chemical properties. Work on this model included specific 
calibration for the Landsat TM sensor. The model treats atmospheric 
optics, water optics, and the wind ruffled air-water interface. A 
solar ephemeral model has also been implemented to provide a capabil- 
ity to simulate the entire sun-sensor geometry. For many of the lakes 
involved in this study absorbing effects of DOC dominate the scatter- 
ing effects of suspended minerals and organic particles. Under these 
conditions subsurface reflectance can be estimated as the ratio of 
backscattered radiation to the total lost by both backscattering (Bb) 
and absorption (a). 

The specific values of a and Bb will depend on the concentrations 
of silt (mineral particles), chlorophyl 1-a pigments (C) , and DOC. The 
absorption and scattering cross sections used in the present study 
were those derived by Bukata [1985] in his detailed optical analysis 
of Lake Ontario waters. These cross sections are shown in Figures 7.1 
and 7.2. 

The specific concentrations of each component were used together 
with these cross sections to estimate the absorption and backscatter- 
ing coefficient. The following equation gives the general subsurface 
reflectance model: 

R < X > - C o< x > • a (X) b i X U) 

where R(A) = Subsurface irradiance reflectance 
C 0 (A) = Constant (typical value = .33) 

Bb(\) = Total backscattering coefficient 
a(\) = Total absorption coefficient 
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Figure 7.1 . Absorption Cross Sections for Chlorophyll-a, DOC, 
Suspended Minerals, and the Absorption Coefficient 
of Pure Water. 





63 


Figure 7.2. Backscatter Cross Sections for Chlorophyll-a, 
Suspended Minerals, and the Backscatter 
Coefficient of Pure Water. 
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This model calculates subsurface reflectances (at the wavelengths 
measured by the MER) given the concentrations of chlorophyll, DOC, and 
suspended solids as shown in the following equation: 

R(X) - C (X) • (Bb„(X) + Bb c (XHC] + Bb^XMSH]) 

° U„(X) + a c (X)*LCJ + a SM (X)>lSMJ ♦ a DOC (X)*lOOCJ + Bb's) 


where R = Subsurface hemispherical reflectance 
SM = suspended solid concentration (mg/1) 

C = chlorophyll concentration (/*g/l) 

DOC = Dissolved organic carbon concentration (mg/1) 


7.2 MODEL CALIBRATION 

Backscattering and absorption values were regressed with the MER- 
1000 estimated subsurface reflectance at each wavelength producing an 
estimate of constant coefficient (C Q ) which is listed in Table 7.1. 
The resulting set of reflectance equations can be used to examine the 
spectral reflectance dependence on DOC and other constituents. The 
mineral particle concentrations were found to be extremely small, on 
the order of 0.1 mg/1. If one assumes a chlorophyll-a concentration 
of 1.0 /ig/1 (a typical value) then the DOC reflectance varies between 
1% and 6% in TM band one as depicted in Figure 7.3. 


7.3 MODEL EXTENSION WITH PR0BAR DATA 

The PR0BAR above-surface reflectance data were collected in August 
1986. These data were converted to subsurface reflectances for over 
one-hundred lakes using a regression procedure (described in Section 
8.5). 

The model developed for the MER subsurface reflectance data was 
tested using the PROBAR-predicted subsurface reflectance data. The 
Marquardt method was used for developing the non-linear model. This 
method is equivalent to performing a series of ridge regressions and 
is most useful when the parameter estimates are highly correlated. 
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Table 7.1. 

Reflectance Model Coefficients 
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Figure 7.3. Reflectance Model for Dissolved Organic Carbon (DOC) 



Since DOC and chlorophyll, (the two model parameters), have a correla- 
tion coefficient of about 0.73, the Marquardt method seemed 
appropriate. 

To estimate how well this model fit the PROBAR predicted subsur- 
face reflectance data, the coefficients produced using these data were 
compared to those produced using the MER data. The results of using 
the non-linear model on data from wavelengths of 443, 470, 520 and 540 
pm are listed in Table 7.2. 

TABLE 7.2. COMPARISON OF PROBAR AND MER MODEL COEFFICIENTS 



PR0BAR 

MER 

443 

470 

520 

550 

.51 

.48 

.42 

.32 

.73 

.68 

.36 

.32 


The model fits the data best in the longer wavelengths. At worst, the 
model coefficients are different by .22, or approximately 30% (for 
X=443 nm) . At best, there is no difference between the coefficients 
(X=550 nm) . 

A comparison of the actual PROBAR predicted subsurface reflectance 
and the model -predicted subsurface reflectance was made to test the 
performance of the reflectance model. The correlation between the 
predicted and actual subsurface reflectance models was quite high, 
ranging from .81 to .89, depending on the wavelength. Model -predicted 
versus PROBAR-predicted subsurface reflectances at 440 nm and 470 nm 
are shown in Figures 7.4 and 7.5, respectively. 

7.4 REFLECTANCE SENSITIVITY TO CHANGES IN WATER CHEMISTRY 

The sensitivity of reflectance to changes in DOC is given by the 
following derivative of the model equation: 
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Fiqure 7.5. Model Prdicted Versus PROBAR Predicted Subsurface 
Reflectance at 470nm. PROBAR Data Collected From 
Algoma and Sudbury Site, August 1986. 
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d R(X) _ a DQC (X).(Bb w (\)+Bb SM (X)»[$M]-»-Bb c (X)»[C]«Co(X) 
d [D°C] (a DOC (X)*[DOC]+a c (X)*[C]+a SM (X)*[SM]+a w (X )) 2 

Figure 7.6 shows the change in reflectance sensitivity for a given DOC 
concentration. The plotted sensitivity values are for the Sudbury 
site, calculated using the above equation and measured values of DOC 
and chlorophyll-a. 

7.5 MODEL-PREDICTED SENSITIVITY OF TM 

The ability to detect a seasonal change using depended on the 
measured TM reflectance changes, and on the sensitivity of reflectance 
to changes in DOC and chlorophyll-a pigment concentration. 

The impact of DOC changes on reflectance can be calculated using 
the sensitivity equation in Section 7.4. The expected TM band one 
signal change per percent subsurface reflectance change was estimated 
previously to be 2.86 counts/percent. If it is assumed that seasonal 
changes in DOC are on the order of 50%, then background levels of two 
to three count changes are projected in the TM response. These pre- 
dictions are summarized as Table 7.3. 


70 


DOC (mg/I) 


Figure 7.6. Sensitivity of Reflectance to Changes in DOC 
Concentration for a Clear Lake Typical of the 
Sudbury Site. 
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8.0 ANALYSIS OF RADIOMETRIC DATA RELATIONSHIPS 

8.1 CHARACTERIZATION OF WATER CHEMISTRY OF STUDY AREA LAKES 

The August 1986 water chemistry data collected in this experiment 
contain twenty-eight in-lake water parameters for 300 lakes across 
Ontario. Pearson correlation coefficients and their significance 
probabilities were produced for a subset of these data set and are 
listed in Table 8.1. There were strong correlations between pH and 
total inflection point alkalinity and aluminum (.88 and -.75, respec- 
tively). The correlation between pH and DOC was found to be much 
lower at 0.61 but which still indicates a significant relationship 
exists. A scatter plot of pH and DOC is shown as Figure 8.1. It is 
evident from these data that the strongest relationship exists for DOC 
values less than 3.0 mg/1. 

8.2 ANALYSIS OF SUBSURFACE IRRADIANCE MEASUREMENTS 

Based upon the reflectance model analysis high correlations were 
expected between lake water chemistry and MER optical measurements. 

The Pearson correlation coefficients and their significant probabili- 
ties are tabulated in Table 8.2 for the August 1986 water chemistry 
data. In general, there is a high correlation between the short wave- 
length reflectances (X < 540 nm) with Secchi depth (SD) , chlorophyl 1 -a 
(CHL0R), DOC, aluminum (AL) , and pH. The high correlations with SD, 
DOC, and AL support the phenomenological relationships between water 
chemistry parameters and optical properties as discussed previously in 
Chapter 2.0. The lower correlations with chlorophyll-a values were 
expected since pigment concentrations measured in many of these lakes 
was so small . 

Mer spectral reflectances were plotted for selected lakes which 
are given as Appendix F. The clear acid lakes were found to have 
spectral reflectances with peaks in the 400-450 nm range and shape 
similar to that obtained for Sunnywater Lake (see Figure 8.2). By 
contrast the high DOC lakes have spectral reflectance curves which 
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Table 8.1 . Pearson Correlation Coefficient for Water Chemistry 

Parameters With Their Significance Probabilities Given 
Directly Below Each Value. 
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TM4M 0.00220 0.15893 -0.07001 -0.20989 0.31086 -0.03466 

0.9855 0.1855 0.5618 0.0789 0.0083 0.7742 




★ ALGOMA 
a SUDBURY 

Figure 8.1. Dissolved Organic Carbon Versus pH Value for 
Water Samples Collected From the Algoma and 
Sudbury Sites, August 1986. 



Table 8.2. Pearson Correlation Coefficient for Water Chemisti 
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R689 -0.76856 0.88722 0.94489 0.92644 -0.66073 -0.76162 

0.0436 0.0077 0.0013 0.0027 0.2001 0.0467 




Figure 8.2. Spectral Reflectance for Sunnywater Lake as Derived 
From MER Data Collected 13 August 1986. 
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Figure 8.3. Spectral Reflectance for Center Lake as Derived 
From MER Data Collected 22 August 1986. 
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have generally lower reflectances values and a spectral peak at 
approximately 550 nm. For these lakes, the high DOC concentrations 
(2.0 - 4.0 mg/1) are consistent with the low reflectance values 
derived for the shorter wavelengths. 

The high-DOC basic lakes have curves shaped more like Center Lake 
(see Figure 8.4). Therefore, the following indicator for character- 
izing acid and basic lakes using these spectral data could be 
calculated: 


j _ Reflectance (500 am) 

" Reflectance (560 /cm) 

This suggested indicator, I, which takes advantage of the difference 
in the shapes of spectral curves, is greater than 1.0 for acidified 
lakes and is less than 1.0 for buffered, high DOC lakes. 

The MER reflectances were also analyzed using the non-linear 
reflectance model described in Section 7.0. The suspended solids were 
assumed to be constant at 0.1 mg/1. The model converged for all the 
MER data collected and the following coefficients C Q (X) are shown in 
Table 8.3. 

TABLE 8.3. COEFFICIENTS FOR SUBSURFACE REFLECTANCE MODEL 
USING MER DATA 

1986 1987 


Wavelength (j»m) 

Aua_ 

May 5 

May 12 

June 13 

June 30 

488 

.524 

.388 

.523 

.779 

.89 

560 

.302 

.338 

.332 

.523 

.667 


The August data were collected under the best conditions, so the coef- 
ficients produced for these data were used as standards to compare the 
other dates. The May 12 data produced coefficients nearly equal to 
those produced using the August data. The June reflectance data do 
not seem to fit the same model suggesting that the water chemistry had 
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Figure 8.4. Comparison of MER and PROBAR Derived 
Spectral Reflectances. 
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changed dramatically and the DOC reflectance model assumptions were no 
longer valid. 

To find out how well the model worked for each date, the correla- 
tions between actual and model -predicted subsurface reflectances were 
calculated. There was no correlation between actual and predicted 
subsurface reflectance for any of the spring data at 560 /<m. For 488 
pm the correlation between actual and predicted reflectances was less 
than .24 for the two June dates. However, the same correlations for 
May 5 and May 12 are .93 and .97, respectively. The reflectance dif- 
ferences between actual and predicted reflectances were less than 
1.15% for these two dates. 

8.3 ANALYSIS OF SURFACE MEASUREMENT DATA 

The PROBAR-derived surface reflectance data were found to be 
highly correlated with the MER subsurface reflectance data as shown by 
the examples in Figure 8.3. To determine if the correlations of water 
chemistry with PROBAR data were similar to those with the MER data, 
another correlation matrix was calculated. Table 8.4 contains PROBAR 
correlations with water chemistry at multiple wavelengths. The cor- 
relations of reflectance with the water chemistry are much lower, but 
still reach -.80, -.68, and -.64 for DOC, pH and chlorophyll. This 
was expected, however, since varying lake surface waves and atmo- 
spheric haze introduced more noise in the signal measured by the 
PROBAR. 

8.4 THE COMPARISON OF SURFACE AND SUBSURFACE MEASUREMENTS 

An experiment was conducted to determine the relationship between 
the surface and the subsurface measurements of lake volume reflec- 
tance. The surface reflectance was measured using the PROBAR spectral 
radiometer mounted in a helicopter and the subsurface reflectance was 
measured using the MER submersible radiometer. Modeling theory pre- 
dicted that the relationship between these two measurements would be 
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Table 8.4. Pearson Correlation Coefficient for W; 
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linear. Therefore, the relationship between PROBAR and MER measure- 
ments could be described using the equation: 

R.j (mer) = m x R^ (PROBAR) + b 

where m = slope 

b = y-axis intercept 
R.j = reflectance at band i 

The results of a linear regression analysis of each spectral band and 
the corresponding statistical significance (a = .05) of each regres- 
sion are found in the following table: 


Wavelength 


Significant 


443 

- .33 

.53 

yes 

470 

.215 

.44 

yes 

490 

.43 

.42 

yes 

520 

.55 

.44 

yes 

550 

1.17 

.19 

no (p = .36) 

580 

1.04 

.14 

no (p = .76) 

610 

-1.11 

1.59 

yes 

640 

- .57 

1.86 

yes 

670 

- .27 

1.49 

yes 

700 

0.0 

1.0 

yes 


The results shown in Section 5.5 support the hypothesis that there 
is a linear relationship between the MER and the PROBAR data for all 
but wavelengths 550 and 580 nm. At most wavelengths, then, subsurface 
lake volume reflectance can be predicted with reasonable accuracy when 
only the PROBAR reflectance data are available. This result is sig- 
nificant since acquiring lake reflectance data with the PROBAR is less 
expensive and quicker that with the in situ measurements. 
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8.5 ANALYSIS OF TM MEASUREMENTS 

The haze normalized TM data for August 1986 show sensitivity to 
lake DOC concentrations as indicated by the data plotted in Figure 
8.5. These results confirm the model predicted sensitivity of the TM 
band on signals to changes in DOC. The model predicted a DOC reflec- 
tance range of about 5% which corresponds to a 14.3 signal count 
spread in the TM band one data. TM data from the Sudbury site are 
consistent with the predicted spread in DN counts. The Algoma data 
appear to lack sensitivity to changes in DOC which is likely due to 
the fact that most lakes in the Algoma region have high values of DOC 
and chlorophyll -a. 

8.6 MULTITEMPORAL RELATIONSHIPS 

Multi temporal analyses could be made for TM and MER data only 
since these were the only measurements made in the spring of 1987. 
PROBAR multitemporal relationships could not be examined since this 
instrument was not available to the project in 1987. 

8.6.1 MER Multitemporal Analysis 

The corrected MER data yielded several multi temporal trends. 

These trends differed depending on the buffering capacity of the lake. 
Acidified lakes, such as Dougherty and Wolf, (TIP < 0), had small 
multi temporal reflectance changes from 500 pm to 600 pm (< .4% reflec- 
tance). All the acid lakes showed large reflectance differences in to 
400 pm to 500 pm region. Each lake showed a decrease in reflectance 
form August, 1986 to May 5, 1987, and then an increase in reflectance 
from May 5, 1987 to May 12, 1987 between 400 pm to 500 pm. These data 
for three lakes' reflectances at 441 pm are tabulated below: 
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Figure 8.5. TM Band 1 Versus Dissolved Organic Carbon Using 

the August 13, 1986 (PI 9, R27) and August 18, 1986 
(P22, R27) Data Sets. 




Name 

8/86 

5/5/87 

5/12/87 

Sunnywater 

6.5 

- 

7.4 

Wolf 

3.2 

2.1 

3.2 

Dougherty 

3.8 

2.3 

2.6 


In contrast, the buffered lakes, (TIP > 0), did not show a large dif- 
ference (> .4% reflectance) in the 400 pm to 500 pm region but did 
show large multi temporal differences from 500 pm to 600 pm. At 
560 pm, the basic lakes increased in reflectance form August 1986 to 
May 5, 1987. No consistent change in reflectance from May 5, 1987 to 
May 12, 1987 was found for the buffered lakes. The reflectance data 
for 560 pm measured from buffered lakes are found below: 

Name 8/86 5/5/87 5/12/87 

Centre 1.8 2.6 1.7 

Whitepine 2 1.2 1.6 2.0 

In conclusion, differences in multitemporal MER reflectance trends 
between buffered and acidified lakes were found. Acidified lakes had 
a decrease in reflectance for the 400 pm to 500 pm region and rela- 
tively no change for the 500 pm to 600 pm region. Buffered lakes had 
an increase in reflectance for the 500 pm to 600 pm region and rela- 
tively no change to the 400 /jm to 500 pm region. 

8.6.2 TM Multitemporal Analysis 

The TM band one seasonal change patterns are similar to those 
indicated for the August 1986 data. The August low DOC lakes were 
found to have larger TM DN values than with the May 12, 1987 and June 
13, 1987 collection dates. These data are shown in Figures 8.6 and 
8.7, respectively. The extracted and atmospherically normalized TM 
data are given as Appendix G. The size of the TM band one count chan- 
ges for Sudbury are substantially larger than predicted. Furthermore, 
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these count differences suggest a two to three percent change in sub- 
surface reflectance, needed a greater DOC change sensitivity than 
predicted in Table 7.3. 

Multi temporal differences were calculated for the following pairs 
of dates: 


August 1986 - May 1987 
August 1986 - June 1987 
June 1987 - May 1987 

These differences were analyzed to determine whether or not they aided 
in identifying acidified and buffered lakes. 

The multi temporal changes in MER-derived subsurface reflectance 
were used to determine the expected changes in TM signal counts for 
bands one and two using the conversion factors given in Section 2.5. 
The expected changes in TM signal counts for band two were all found 
to be within the noise level for band two data. Therefore, the band 
two multitemporal differences were insignificant for the 1986-1987 
scene pair. Furthermore, approximately 90% of the TM band 1 dif- 
ferences were also in the noise level. As a result, obvious multi- 
temporal differences using TM data were not found. 

The majority of the multitemporal analyses were based on the 
August - May scene pair. When all of the lakes (n = 41) are analyzed 
for significant differences (a = .1) between August and May reflec- 
tance changes, no difference is found between acidified and buffered 
lakes. 

Another test was made to determine if the August - May TM band one 
differences were a function of DOC, TIP, chlorophyll and/or aluminum 
levels measured in 1986. A multivariate regression analysis was done 
and all of the parameter estimates were insignificant (a = .5). These 
results lead to conclusion that the TM band one differences were not a 
function of water chemistry. 


88 



Data were also extracted from a May 22, 1985 scene for the Sudbury 
site and the differences in TM band 1 were formed with the August 13, 
1986 scene. These differences were then compared to DOC values col- 
lected in August 1986. The results are shown in Figure 8.8 and indi- 
cated a possibly strong seasonal relationship to DOC concentrations 
and especially for those lakes with less than 3.0 mg/1 DOC. 

8.7 ANALYSIS OF TRANSMISSOMETER ATTENUATION DATA 

A study was conducted which examined the relationship between the 
water attenuation coefficient (a) at 600 /*m and the suspended solid 
concentrations in eight lakes. The attenuation coefficient correlates 
positively with the suspended solids ( p = .845). These data are shown 
in Figure 8.9 and shows and a linear relationship between the attenua- 
tion coefficient and the suspended solids. This further supported the 
possibility of suspended solid concentrations affecting the accuracy 
of the subsurface reflectance model since the suspended solids concen- 
tration correlates with the attenuation coefficient and the attenua- 
tion coefficient affects lake volume reflectance. 
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Figure 8.8. TM Band 1 Multitemporal (August 13, 1986 and May 22, 1985) 
Differences Versus DOC Concentration Sudbury Field Site 
August 1986 Water Chemistry Data. 
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Figure 8.9. Beam Attenuation Coefficient Versus Suspended Solids 
Concentration 1 987 Spring/Summer Data. 
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9.0 ANALYSIS OF ECO-PHYSICAL CLUSTERS 

9.1 RELATIONSHIP OF WATER CHEMISTRY WITH ECO-PHYSICAL CLUSTERS 

Since the eco-physical clusters represented unique acid- 
sensitivities across the Ontario test sites, it is reasonable to 
expect to find significantly different water-quality parameters for 
lakes that occurred in different clusters. An analysis was performed 
to test the hypothesis that the mean water-quality parameters were 
different at the 5% significance level between clusters. The follow- 
ing water variables were analyzed: dissolved organic carbon (DOC), 

Secchi depth, sulfate concentration, aluminum ion concentration, pH 
and total chlorophyl 1-a concentration. This analysis is summarized in 
Table 9.1. 


TABLE 9.1. RESULTS FOR TUKEY'S STUDENTIZED RANGE TEST FOR 
SIGNIFICANTLY DIFFERENT MEAN WATER-QUALITY 
PARAMETERS 


Chloroc 


Secchi Depth Sulfate Aluminum 


1-7 2-5 1-7 

1-9 2-9 2-7 

1-8 2-7 2-9 

1-5 4-5 4-7 

4-7 5-7 

7-5 
4-9 
3-9 
3-7 

1- 7 

2 - 8 

3- 5 

4- 8 
7-8 


7-1 

7-5 

1-2 

8-2 

7-2 

3-7 

1-3 

8-3 

7-3 

4-5 

1-4 

8-5 

7-4 

5-9 

1-5 

8-7 

8-9 

7-8 

1-6 

3-9 

1-2 


1-7 

2-7 

7-8 


1-8 

4-2 

7-9 


1-9 

4-3 

9-1 


1-10 

4-5 

9-3 


6-2 

4-7 

9-4 


6-3 

4-9 

5-1 


6-5 

3-2 

5-4 


6-7 

3-7 

3-1 


6-9 

5-7 

1-4 


10-2 

9-7 

1-6 


10-3 


5-7 


10-5 




10-7 




10-9 



PRECEDING PAGE BLANK NOT FILMED 

93 


past 1 ± WTiwiONMt? turn 




Clustering was especially successful in separating different levels of 
lake pH. For significantly different clusters, the most acid- 
sensitive clusters (5,7,9) had lake waters with lower DOC and pH 
values and a higher sulfate concentration than those with less sensi- 
tivity (1,2, 3, 4). Thus, the clustering analysis appears to have pro- 
duced significant eco-physical clusters that contain lakes that also 
have some significantly different water-quality parameters. Of the 
three most significantly different water-quality parameters (DOC, pH 
and sulfate), changes in DOC provide the basis for remote sensing 
monitoring and identification. 

9.2 RELATIONSHIP BETWEEN TM SIGNALS AND ECO-PHYSICAL CLUSTERS 

An analysis was performed to determine if significant differences 
existed among the eco-physical clusters based on the TM signals of 
lakes within the clusters. For the August 1986 Algoma and Sudbury 
data sets, two groupings were identified. Group A (lakes in clusters 
5, 7 and 9) had mean signals (73.5 to 75.9) that were significantly 
different (a= .05) than signals (64.8 to 67.5) from lakes in group B 
(clusters 2, 4 and 8). The results are shown in Table 9.2. 

TABLE 9.2. TM RELATIONSHIPS TO ECO-PHYSICAL SENSITIVITY — 
AUGUST TM 1 DATA 


Group 

Mean TM 1 

Cluster 

Siqnificantly Different at 5% 

A 

75.86 

7 

00 

** 

C\l 

A 

74.14 

9 

2, 4, 8 

A 

73.47 

5 

2, 4, 8 

B 

67.50 

8 

5, 7, 9 

B 

66.37 

2 

5, 7, 9 

B 

64.77 

4 

5, 7, 9 


The mean eco-physical sensitivity of group A clusters was 7.44 and 
group B mean sensitivity was 5.85. The largest signals measured were 
from cluster 7 with a mean sensitivity index of 7.07 and the smallest 
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from cluster 4 with a sensitivity index of 6.07. The primary dif- 
ference in these two eco-physical clusters was the soil type and soil 
depth over the underlying bedrock. Cluster 7 had shallow (i.e. less 
than one meter) sandy soils while cluster 4 had soils of mixed types 
(sand, clay, loam) that had depths greater than one meter. 

9.3 RELATIONSHIP BETWEEN TM MULTITEMPORAL DIFFERENCES AND ECO- 
PHYSICAL CLUSTERS 

Examination of the August 1986 - May 1985 difference signals for 
TM band one produced similar results which are shown in Table 9.3. 
Group A and group B contained the same clusters as in Section 9.2 and 
the largest and smallest mean differences were found in clusters 7 and 
4, respectively. 


TABLE 9.3 


TM RELATIONSHIPS TO ECO-PHYSICAL SENSITIVITY 
ANALYSIS OF VARIANCE OF AUGUST-MAY DIFFERENCE 


Group 

Mean Diff 

Cluster 

Sianificantlv Diff. at 5% 

A 

4.94 

7 

2, 4, 8 

A 

2.91 

9 

2, 4, 8 

A 

2.68 

5 

00 

*1- 

CM 

B 

0.61 

8 

o> 

in 

B 

-0.96 

2 

5, 7, 9 

B 

-2.45 

4 

5, 7, 9 


9.4 Analysis of TM Signal Changes Due to Acid Deposition Changes 

This analysis examined the relationship between the August - May 
signal differences from polygons that have similar eco-physical prop- 
erties with the exception of sulfate deposition. For this case, lakes 
were selected from polygons with sandy soils over granitic rock types 
and the sulfate deposition was 1.5 or 2.5 g/m2/yr. The TM band one 
signals were found to be significantly different (at 5-tf level) based 
upon deposition level alone. This preliminary analysis suggests that 
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TM signal seasonal changes may be dependent upon changes in acid 
deposition. 

9.5 ANALYSIS OF DOC REFLECTANCE SENSITIVITY 

In addition to seasonal analyses, the spatial aspects of DOC 
reflectance sensitivity were investigated. Measured water-quality 
parameters were used together with the equation given Section 7.4 
to calculate a lake value of reflectance sensitivity based to change 
in DOC concentration. The larger the derivative of reflectance with 
respect to DOC the more sensitive lake reflectance is to changes in 
DOC. As shown in Table 7.3, a reflectance sensitivity value of 4.0 
corresponds to an expected count change in TM band 1 of at least two 
counts. The lake DOC sensitivity values were analyzed with the eco- 
physical clusters and the mean sensitivity was determined for each 
cluster as shown in Figure 9.1. These results indicated that clusters 
5, 7, 8 and 9 have lakes most sensitive to DOC changes. These clus- 
ters also have the higher stratification sensitivity index values. 

This preliminary analysis shows that TM band one seasonal differ- 
ence signals will differentiate acid-sensitive from acid-insensitive 
areas. 
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Figure 9.1. Mean DOC Induced Reflectance Sensitivity for 
Each Ecophysical Strata Estimates Based Upon 
August 1986 Water Chemistry Measurements. 
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10.0 CONCLUSIONS AND RECOMMENDATIONS 

10.1 GENERAL CONCLUSION 

Modeling, field measurements, and TM observations suggest that TM 
is useful to identify acid sensitive lakes and to monitor water qual- 
ity changes associated with lake acidification. 

10.2 SPECIFIC CONCLUSIONS 

1. Modeling surface and subsurface reflectance measurements have 
shown the important relationships between lake optical proper- 
ties and water chemistry. 

a. A simple DOC reflectance model accounts for observed sub- 
surface hemispherical reflectance and also for the compan- 
ion airborne (PROBAR) reflectance measurements. 

b. We found that clear acid sensitive lakes can be distin- 

f 

guished from the colored high DOC lakes using PROBAR data. 

The colored lakes tend to have greater buffering capacity 
than clear lakes in acid sensitive areas. 

c. The blue-green reflectance of clear lakes is highly sensi- 
tive to the presence of DOC. Modeling predicts a one 
percent change in subsurface reflectance for an expected 
seasonal fluctuation of about 50% in DOC concentration. 

2. Modeling has shown that TM is sufficiently sensitive to moni- 
tor expected lake reflectance associated with acid deposition 
and acidification. 

a. An historical TM seasonal pair (August 1986 - May 1985) 
for the same Sudbury Lakes in a normal snowfall year sup- 
ports our expectations but lacks the chemistry and in situ 
optical data needed for hypothesis validation. 

b. The expected seasonal changes (August 1986 to May/June 
1987) in water chemistry did not occur nor did we observe 

PRECEDING page BLANK NOT F(LMED 


PA 6 I ‘Vj J HItffltQNAmi WAHK 


a significant change in TM response. This lack of change 
may be due to the unusually small snow pack and spring 
runoff. While these TM data and water chemistry are con- 
sistent with our hypothesis, they do not validate it. 

c. In areas of high acid deposition Landsat TM DN values were 
found to separate high DOC lakes (moderate acidity) found 
to separate high DOC lakes (moderate acidity) from low DOC 
lakes (high acidity). The expected TM band one signal 
change per percent subsurface reflectance change was esti- 
mated to be 2.86 counts/percent. In clear acid lakes 
seasonal change of two or three counts are expected from 
DOc fluctuations. 

Stratification of eco-physical properties provides a way to 

locate areas which are sensitive to acid deposition. 

a. When stratification of eco-physical properties was applied 
to our study sites, we could identify acid sensitive areas 
and use TM to pick lakes which are likely to be sensitive 
to acid deposition. 

b. Clustering of eco-physical strata suggests that areas with 
shallow sandy soils over slow weathering granitic bedrock 
types are most sensitive to acid deposition and lakes 
located within these areas will have lower concentrations 
of DOC and lower pH values than for other soil and bedrock 
types . 

c. TM band one lake response was found to be related to eco- 
physical sensitivity. The (August 1986 - May 1985) TM 
seasonal pair produced signal differences in eco-physi- 
cally sensitive strata (1-6 DN) but not so in non-sensi- 
tive strata (-2 to 0 DN) . 

d. Nearly identical and sensitive eco-physical strata with 
different sulfate deposition rates were found to have 
different TM lake signal response. 
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10.3 RECOMMENDATIONS 

While studies thus far are consistent with our seasonal change 
hypothesis they do not confirm its validity. Further study is needed 
to provide confirmation to the above results. 

1. Collect lake chemistry and TM data in years of typical snow- 
fall to demonstrate the capability of using TM data to monitor 
acidification under a wider range of environmental conditions 
(i.e., normal snowfall years). 

2. Develop a TM based capability for assessing effects of acid 
deposition on terrestrial vegetation. Apply the vegetation 
monitoring technique and compare with lake monitoring tech- 
nique. 
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APPENDIX A 

ECO-PHYSICAL CLUSTER ANALYSIS 

The maximum likelihood method was used to produce 10 clusters of 
polygons based on the sensitivity values for percent cover, vegetation 
type, soil depth, soil texture, bedrock type, relief and sulfate depo- 
sition. The data are sorted by cluster. Descriptions for each poly- 
gon are in the printout. The "cluster" data are either missing data 
or have vegetation types which were not used in the data. 
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ERIM 


APPENDIX B 

PR08AR REFLECTANCE DATA 

Table B.l. Corrected PROBAR reflectances above the water surface and 
water chemistry data. 

Table B.2. PROBAR subsurface predicted reflectances. 
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PREDICTED SUBSURFACE REFLECTANCE AT DEPTH = 2 METERS 
ALL PROBAR LAKES 
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APPENDIX C 

SUMMARY STATISTICS FOR THE ECO-PHYSICAL POLYGON CLUSTER ANALYSIS 


The following table shows the computer mean and standard deviation 
estimates for the set of Eco-physical polygons within each cluster. 
Estimates are computed for the total sensitivity rating (STRATRAT) , 
vegetation sensitivity (VEGVAL) , bedrock and soil sensitivity 
(SENSVAL) , relief sensitivity (RELVAL), and sulfate deposition sensi- 
tivity (S04VAL) . 
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TABLE C-l 


SUMMARY STATISTICS ON EACH CLUSTER 
MAXIMUM LIKELIHOOD CLUSTER ANALYSIS 


VARIABLE 

MEAN 

STANDARD DEVIATION 

Cluster= 

STRATRAT 

SENSVAL 

3.13 

0.71 

VEGVAL 

1.42 

0.69 

RELVAL 

0.38 

0.12 

S04VAL 

1.30 

0.24 

Cluster=l 

STRATRAT 

5.66 

0.07 

SENSVAL 

7.04 

0.66 

VEGVAL 

4.68 

0.75 

RELVAL 

5.57 

0.23 

S04VAL 

4.40 

0.33 

Cluster=2 

STRATRAT 

6.36 

0.06 

SENSVAL 

8.05 

0.42 

VEGVAL 

4.65 

0.66 

RELVAL 

5.78 

0.20 

S04VAL 

5.82 

0.38 

Cluster=3 

STRATRAT 

6.74 

0.09 

SENSVAL 

8.16 

0.32 

VEGVAL 

5.83 

0.64 

RELVAL 

5.28 

0.22 

S04VAL 

6.00 

0.36 

Cl uster=4 

STRATRAT 

6.02 

0.07 

SENSVAL 

7.67 

0.52 

VEGVAL 

4.63 

0.67 

RELVAL 

5.25 

0.21 

S04VAL 

5.18 

0.42 
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Cluster=5 


STRATRAT 7.41 
SENSVAL 8.47 
VEGVAL 7.13 
RELVAL 5.62 
S04VAL 6.58 


0.06 

0.20 

0.47 

0.19 

0.29 


Cluster=6 


STRATRAT 3.55 
SENSVAL 3.28 
VEGVAL 2.08 
RELVAL 5.57 
S04VAL 5.27 


0.29 

0.92 

0.14 

0.17 

0.68 


Cluster=7 


STRATRAT 7.07 0.05 
SENSVAL 8.50 0.21 
VEGVAL 6.37 0.48 
RELVAL 5.36 0.22 
S04VAL 6.10 0.32 


Cluster=8 


STRATRAT 

5.14 

0.20 

SENSVAL 

5.96 

0.57 

VEGVAL 

4.71 

0.59 

RELVAL 

5.46 

0.22 

S04VAL 

3.97 

0.33 


Cluster=9 


STRATRAT 

7.83 

0.20 

SENSVAL 

8.72 

0.22 

VEGVAL 

8.53 

0.49 

RELVAL 

5.20 

0.22 

S04VAL 

6.30 

0.29 


C1uster=10 


STRATRAT 

4.34 

0.22 

SENSVAL 

5.22 

0.29 

VEGVAL 

3.82 

0.40 

RELVAL 

5.00 

0.22 

S04VAL 

3.05 

0.21 
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APPENDIX D 

WATER CHEMISTRY DATA 


Table D.l. August 1986 WQ Data Collected from the Algoma and Sudbury 
sites 

Table D.2. May - June 1987 WQ Data Collected from selected lakes in 
the Sudbury site 

Figure D.l MER and PROBAR Sampling Stations for the Algoma Site 

Figure D.2 MER and PROBAR Sampling Stations for the Sudbury Site 

Maps shown in Figures D.l (80798) and D.2 (80799) were compiled by 
J. Fortescue and D. Stahl of the Mines and Minerals Division, Ontario 
Geological Survey, 1987. 
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AUGUST 1980 WATER CHEMISTRY 
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SPRING 1987 WATER CHEMISTRY 
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Figure D.l. MER and PROBAR Sampling Stations 
for the Algoma Site 
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APPENDIX E 

TRANSMISSOMETER DATA DERIVED TRANSMISSION AND ATTENUATION COEFFICIENTS 
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SEA TECH TRANSMISSOMETER 
SUMMER AND SPRING DATA 
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#ERIM 


APPENDIX F 

MER-SUBSURFACE SPECTRAL RADIOMETER MULTITEMPORAL LAKE REFLECTANCES 


Figure F.l 

Smoothwater Lake 

Figure F.2 

Whitepine #1 Lake 

Figure F.3 

Sunnywater Lake 

Figure F.4 

Wolf Lake 

Figure F.5 

North Yorkston Lake 

Figure F.6 

Whitepine #2 Lake 

Figure F.7 

Dougherty Lake 

Figure F.8 

Centre Lake 
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Figure F.1 Smoothwater Lake 




Smoothwatar Laka Mar Data at 2 Matars 
Multi tampora I 
Mar Raflactanca 
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Figure F.2 Whitepine #1 Lake 
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Figure F.3 Sunnywater Lake 
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Wolf Laka Mar Data at 2 Matars 
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Figure F,5 North Yorkston Lake 


North_Yorston Lake Mar Data at 2 Matars 
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Whitapina_#2 Laka Mar Data at 2 Matars 
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Figure F.8 Centre Lake 




Caotra Lika Mar Data at 2 Matara 
Multi tampora I 
Mar Raflactanca 
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APPENDIX G 

LAKE EXTRACTED TM SIGNAL VALUES AND ATMOSPHERIC CORRECTED VALUES 

Thematic mapper (TM) signal digital count value for lake extracted 
samples are listed in the following tables. Also listed are the 
standard deviation estimates for each sample. 

Table G.l. August 13, 1986 (P19.R27) 

Table G.2. August 18, 1986 (P22.R27) 

Table G.3. May 12, 1987 (P19.R22) 

Table G.4. June 13, 1987 (P19.R27) 

Atmospherically normalized value are listed in the following 
tables. 

Table G.5. August 13, 1986 (P19, R27) 

Table G.6. August 18, 1986 (P22.R27) 

Table G.7. May 12, 1987 (P19.R27) 

Table G.8. June 13, 1987 (P19.R27) 
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AUGUST 18, 1988 
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SUDBURY QUAD 3 
JUNE 13, 1987 

RAW TM SIGNALS AND STANDARD DEVIATIONS 
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Table G.5 

SUDBURY QUAD 3 
AUGUST 13, 1986 

CORRECTED TM SIGNALS AND STANDARD DEVIATIONS 
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Table G.6 

ALGOUA QUAD 4 
AUGUST 18, 1988 

CORRECTED TM SIGNALS AND STANDARD DEVIATIONS 
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